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Abstract 

We study a lattice fermion model for superconductivity in the presence of an 
antiferromagnetic background, described as a fixed external staggered mag- 
netic field. We discuss the possibility that our model provides an effective 
description of coexistence of antiferromagnetic correlations and superconduc- 
tivity, and its possible application to high temperature superconductivity. We 
also argue that, under certain conditions, this model describes a variant of the 
periodic Anderson model for heavy fermions. Using a path integral formula- 
tion we construct mean field equations, which we study in some detail. We 
evaluate the superconducting critical temperature and show that it is strongly 
enhanced by antiferromagnetic order. We also evaluate the superconducting 
gap, the superconducting density of states, and the tunneling conductivity, 
and show that the most stable channel usually has a (4.2 -wave gap. 
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I. INTRODUCTION 



Motivated by heavy-fermion- and high-temperature superconductors, considerable the- 
oretical effort has gone into explaining the occurrence of unconventional superconductivity in 
strongly correlated electron systems. A basic model of interest here is the one-band Hubbard 
model which describes lattice fermions interacting with a strong on-site Coulomb repulsion 
0. This model is a widely studied candidate for explaining the unusual properties of the 
high temperature superconductors. Despite massive efforts this model is not yet well under- 
stood due to the technical difficulties involved. It is thus desirable to devise approximate 
simplified models focusing on certain aspects of the physics, that are easier to analyze than 
the full model, and whose properties can be compared with experiments and simulations. In 
this paper we study such a model focusing on the possible influence of antiferromagnetism 
on superconductivity. 

There are several reasons for studying the coupling between antiferromagnetic (AF) 
correlations and superconductivity (SC): (i) The half filled Hubbard model is known to 
be an insulating antiferro magnet, and quantum Monte Carlo results |l| and results from 
Hartree-Fock theory indicate that local AF correlations persist up to quite high doping 
levels (at least for parameters adequate for high temperature superconductivity), (ii) The 
parent compounds of the high-temperature superconductors (HTSC) are Mott insulators 
with antiferromagnetic long range order. Neutron scattering experiments || show that 
the antiferromagnetic correlation length gets reduced upon doping. NMR measurements, 
however, suggest that strong short range AF correlations remain important all the way into 
the SC phase (see Ref. Q for a recent discussion on this point), (iii) As we will further 
discuss below, AF correlations may also play a role in the periodic Anderson model for 
heavy fermions ||. (iv) Finally there are other (low-T c ) superconductors which exhibit 
simultaneous SC and AF order H0. 

Several different approaches to consider the effects of AF correlations on HTSC have 
been explored in the literature. One frequently studied model is based on fermion pairing 
due to exchange of AF spin fluctuations ||. Another is the spin bag mechanism which 
is based on the idea that local suppression of the AF gap leads to an effective attraction 
between certain quasiparticles, and thus to superconductivity ||. Recently Dagotto and 
coworkers proposed and studied weak coupling models of holes in an antiferromagnet, using 



a numerically determined dispersion relation ||10|| . Most of these models lead to d x 2_ y 2 
pairing, which indeed seems to be observed in experiments on HTSC [11 . 

The model studied in this paper describes lattice fermions in an external antiferromag- 
netic field coupled to the fermion spin. We also include various attractive, instantaneous 
fermion-fermion couplings (see Sec. |11 A[ ). Although we use a field with AF long-range order, 
we argue below that, for the properties we calculate, this model gives a good approximation 
also when the AF order parameter is slowly varying in the sense further discussed below. 
Our model illustrates in a simple way how correlations can narrow the quasiparticle bands 



which increases the density of states (DOS) and hence increases T c [0,0. This is an ex- 



ample of the van Hove scenario which has been frequently discussed in the context of HTSC 
|14] , |10|| . In our model the enhancement of the van Hove singularity by antiferromagnetism 
increases T c , such that the right order of magnitude for HTSC can be obtained by reasonably 
small couplings. Recent numerical results suggest that in the parameter range of interest for 
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HTSC, the Hubbard model has a DOS peak, and the chemical potential intersects this peak 
for certain finite doping [fL5|| . While it is clear that the simple band structure obtained in 
our model is not going to be realistic in detail, our model attempts to (i) represent the states 
at this peak by the AF bands, and (ii) to model the resulting superconducting properties. 
Since our model is simple, our investigation can be done in some detail which would be 



difficult for more realistic models. Our approach is similar to Ref. |10| in that we attempt 
to model quasiparticle bands renormalized by strong correlations and then investigate their 
superconducting properties, but otherwise the models are quite different. 

The model we study naturally leads to a superconducting gap A with a staggered real 
space structure (i.e., with a two component order parameter). A superconducting order 
parameter with such a structure has been studied previously on a phenomenological level in 



the context of heavy fermions |16j . Our results show how such an order parameter can arise 
in microscopic models. 

The first purpose of this paper is to provide motivation for our model and develop efficient 
tools to study it. In any case, our model should describe systems where antiferromagnetic 
order and superconductivity coexist 0. However, we also argue that it can be regarded as 
a useful toy model for studying the effects of band renormalization by strong correlations 
on superconductivity in HTSC or heavy fermion systems. The idea is that by using a 
Hubbard-Stratonovich transformation, an interacting fermion systems can be represented 
as non-interacting fermions coupled to a dynamical boson field. A simple approximation is 
to replace the dynamical boson field by a fixed boson field, which is determined by Hartree- 
Fock (HF) equations. A general feature of the doped Hubbard model is that this HF field 
has a nontrivial spatial dependence 0. This leads to a variety of complicated bands. The 
model studied in this paper represents the simplest nontrivial case where this occurs. To 
study this model, we use an efficient path integral formalism which should be useful also 
in more complicated cases, and we derive mean field equations of BCS type. To further 
investigate the relevance of the description provided by our model a more elaborate analysis 
(e.g. fluctuation calculations) is required, but this is beyond the scope of this paper. 

The second purpose of the paper is to calculate some physical properties of the model. 
Using parameters motivated by HTSC, we present a systematic analysis of the different 
possible SC channels. We study the effect of on-site, nearest neighbor (nn) and next-nearest 
neighbor (nnn) interactions. We find that the dominating channels are neither translation- 
nor spin rotation invariant. In the presence of AF correlations, we find that on-site attraction 
never leads to significant pairing and nn attraction is most efficient and usually produces a 
dj,2_ y 2-wave gap, A oc cos(/ci) — cos^)- For dominant nnn-attraction we get p-wave pairing, 
A oc sin(ki + /C2) + sin^! — fc 2 )- d x 2_ y 2-wave superconductivity has been found in numerous 
other calculations |L1], but our result shows that it arises already in a particularly simple 
model. We also derive and solve the gap equations for nn coupling in the full temperature 
range below T c . Our formalism is convenient to evaluate SC properties. As an example, 
we evaluate the SC density of states and the resulting tunneling conductance as a function 
of temperature. We note that some features of our results appear qualitatively similar to 
HTSC, e.g. the doping dependence of T c and the tunneling conductance. 

The plan of this paper is as follows. In Section [IT] we give a detailed description of the 
model, and then discuss the possibility that it may provide a simple description of certain 
aspects of the physics of the Hubbard model with an additional weak pairing interaction. 
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We also relate our model to the periodic Anderson model for heavy fermions || Section P C . 
Our analysis of the model starts in Section 111 , where we summarize the path integral formal- 
ism. The mean field equations are obtained in Section ffVl We first (Sec. [IV Jj| ) obtain the 
T c -equations for all possible superconducting channels (singlet, and triplet, translational- 
invariant and staggered, s-, p- and (i-wave) which are needed in our systematic stability 
analysis, and then (Sec. [IV C|) derive the gap equations below T c for the dominating chan- 
nels. Section [V] contains our numerical results of our stability analysis, the temperature 
dependence of the SC gap, and the density of states and the tunneling conductance as a 
function of T < T c . We end with conclusions in Section [VI]. Some technical details are in 
two appendices. 



II. THE MODEL 

In this section we present our model (A). We then discuss the relation of our model to 
the doped Hubbard model (B) and the periodic Anderson model (C). 



A. Model description 

Our model Hamiltonian describes lattice fermions with field operators ai + ) (x) in an ex- 
ternal staggered magnetic field B(x) and an additional weak attractive BCS-like interaction, 

#eff = #hop + J2 S W ' B ( X ) e * Q ' X + H ™ W 
x 

where B = (Bi, B 2 ,B 3 ) is the staggered magnetic field representing the AF order, 

H hop = - ta+{x)a a {y) (2) 

is the usual hopping Hamiltonian (i.e. the sum over all nearest-neighboring sites (x, y)), 
S = a + (x)<xa(x) is 2 times the electron spin operator (cr are the Pauli spin matrices), and 

H int = ^^2n(x)V(x-y)n( y y), n = S ^a^a (7 (3a) 

is a (weak) attractive charge-charge interaction. We consider the model on a rf-dimensional 
cubic lattice Z d (i.e. spatial vectors are x = (xi, . . . , Xd) with Xi integers), and Q = (ir, . . . , 7r) 
is the usual AF vector. We use a pairing potential with the Fourier transform 

V(k) = -90 - 9nn ^ 2 COS ( k i) ~ fi'nnn ^ 2 COS (^) COs(%) (3b) 

i i^j 

where k = (ki, . . . , kd) with —n < k{ < 7T and gi > 0. The first, second and third terms here 
describe on-site, nearest neighbor (nn) and next-to-nearest neighbor (nnn) instantaneous 
interactions. We take the corresponding couplings gi as input parameters of our model. 
A useful way to think about this potential is as coming from a Taylor expansion of some 
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general potential V(x — y), where only the lowest order terms are kept since they are the 
only ones expected to be important for SC. Note that the normalization in Eq. (|3b|) is such 
that V(x — y) is — g , —g nn and — g nnn if x and y are equal, nn, and nnn, respectively. We 
do not fix dimension in our formal manipulations, but in our numerical calculations we take 
d = 2. 

We will analyze the model setting B(x) independent of x. This case is of interest for 
materials where AF long range order coexisting with SC has been observed ||7|||. However, 
we expect that our results are also adequate for the case of varying staggered magnetic fields 
B(x), as long as this variation is slow. Arguments for this will be given below. In particular, 
even though the second term in Eq. ([!]) explicitly breaks spin rotation invariance, this model 
can nevertheless describe also situations where this symmetry is unbroken: one can think 
of Eq. (|I|) as coming from an adiabatic approximation (similar to the Born-Oppenheimer 
approximation) of a strong-coupling model in which the dynamical AF fluctuations are 
frozen, and spin rotation invariance will be restored after averaging over certain B(x)- 
configurations. 

B. On the relation to the doped Hubbard model 

Here we describe how our model can be related to the doped Hubbard model with an 
additional attractive interaction H[ nt as in Eq. (|3a|) . The idea is to replace the strong 
Coulomb interaction by coupling the fermions to dynamical bosons. In the path integral 
formalism this can be done by a Hubbard-Stratonovich transformation which introduces 
dynamical boson fields [|T7| 

o (r, x) = r(r, x), 0(r, x) = B(r, x) e iQ x (4) 

depending on imaginary time r. These fields have a simple physical interpretation: s — |B| 
and e = B/|B| represent the magnitude and direction of the AF ordered fermion spins and 
r the fermion charge. 

A standard approximation then is to replace this dynamical boson fields by a fixed and 
static boson field configuration. This formally amounts to a saddle point evaluation of 
the exact boson path integral. The stationary points in this path integral are determined 
by HF equations, and this procedure corresponds to mean field theory. For the doped 
Hubbard model this approximation is already quite involved: Numerical results show that 
the HF solutions have a complicated spatial dependence and, depending on parameters, 
describe magnetic domain walls, magnetic vortices etc. 0. However, in the parameter 
regime of interest all these solutions exhibit antiferromagnetic correlations. The non-trivial 
assumptions required to establish a connection of our model to the Hubbard model are 
(i) HF theory is adequate to determine interacting fermions bands, (ii) it is mainly these 
antiferromagnetic correlations which determine these bands close to the Fermi surface. 

Let us be more specific on point (ii): For thin domain walls or localized vortices (e.g.), 
r(x) = r and B(x) = B are constant in large regions of configuration space. Then there 
are fermion states which, for most x, coincide with states describing a system with constant 
r(x) = r and B(x) = B. Theses states can be viewed as scattering solutions, e.g., in a 
domain wall background. Our model tries to capture the effect of AF correlations on such 
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delocalized states. Our assumption is that such states exist and that for suitable doping, they 
are at the Fermi energy. Note that the filling p used in our model should not be interpreted 
as filling in the Hubbard model: p — 1 is the doping of the AF band, which is smaller than 
the total doping if part of the charge carriers are bound, e.g., in domain walls. A check of 
the properties we assumed requires a rather detailed analysis of the Hubbard model itself 
and is beyond the scope of this paper. These properties seem however to have some support 
from recent numerical simulation results for the Hubbard model, which indicate that for 
parameters adequate for HTSC, p indeed intersects a peak in the DOS at a certain doping 
level |T§. 

We analyze the model setting B(x) independent of x = (r, x), but we expect that our 
results are a good approximation also if B(x) varies slowly or only in small fractions of 
spacetime. Then the 2-point Green functions of our model can be approximated by 

G(x, y) « G N (x -y,y) + G A (x - y, y) e' Qy (5) 

where Gn,a depends on the second argument y only via e(y) ■ er and s(y). These Green 
functions can be obtained from the ones for B(y) = se = const in leading order of a 
gradient expansion. In particular, thermodynamic properties only depend on s(y) in this 
approximation. If s(y) ~ const, we thus expect that the thermodynamic properties of our 
model should not depend much on the AF correlation length (assuming the latter is large 
enough). If s(y) also varies in small fractions of spacetime, the thermodynamic properties 
are appropriate averages of the 'local' properties which can be deduced from the results in 
this paper. 

C. On the relation to the periodic Anderson model 

Here we present an argument that our model, with constant B, can be obtained from 
the periodic Anderson model in a certain parameter regime. The physical picture is as 
follows: the periodic Anderson model is a 3D lattice fermion model with two kinds of spin-| 
fermions: /-electrons with a strong Hubbard repulsion are coupled weakly to non-interacting 
a-electrons 0. For large on-site repulsion, the /-electrons can be regarded as half filled and 
antiferromagnetically ordered. The lower magnetic /-band is then far underneath the Fermi 
surface and the upper one far above. Thus the /-electrons are in a half-filled band and thus 
act as a commensurate [e(x) = const] AF background for the weakly coupled a-electrons, 
and our model provides an effective weak-coupling description of these a-electrons. 

We now turn to a more detailed argument. We start from the following periodic Anderson 
Hamiltonian 

H = H hop + H f + V fa + H mt (6a) 
where Hh op and H int are as above, and 

H f = E("W( X ) + ^/-tWnwW) (6b) 

X 

describes another species of electrons denoted as / which are localized and have a strong 
on-site charge-charge repulsion U » ej » t (nj ]Cr (x) = /+ (x)/ CT (x), and nj = %,t+ w /,j,)- 
The interaction of the a and / electrons is local, 
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Vfa = A£(a+(x)/„(x) + /+(x)a ff (x)) , (6c) 



and i? int is as in Eq. ( |3aj ) i.e., we assume the BCS attraction among the a electrons only. 
H a + Hf + V/ a is the Anderson Hamiltonian: Hf is the usual strongly repulsive Hubbard 
Hamiltonian (with hopping constant tf — ► 0) which in the half- filled case (one /-electron per 
site) describes antiferromagnetically ordered electrons. The a-electrons are assumed in the 
conduction band weakly coupled to the /-electrons, and since U >> e/ >> A 2 and e/ >> t, 
one still expects the / electrons to form a half-filled antiferromagnet. 

The simple physical picture leading to the effective Hamiltonian H e R is that the coupling 
between the / and a electrons and the antiferromagnetic ordering of the former produces an 
effective staggered magnetic field that couples to the spin of the a electrons. The staggered 
magnetic field B is determined by the interband coupling (see below). This effect is inde- 
pendent of the filling in the a band, because the antiferromagnetism comes from the / band, 
which is half-filled. This makes the situation simpler than in the non-half-filled repulsive 
Hubbard model, where commensurate order cannot be assumed. 

In the remainder of this section, we derive if e ff from H by an argument patterned after 
Anderson's derivation of the Heisenberg antiferromagnet from the half-filled Hubbard model. 

We do a second-order perturbation expansion in A to show that the coupling between 
the a and the / electrons produces the staggered field term of the effective Hamiltonian. For 
A = 0, the ground state wave function factors into the product of that of the ground state 
for the a and the /. Since the / system is half-filled because U is large, the ground state 
wave function for the / electrons is 

|Neel) = n/ CT + (x )W|0), a(x)= e ^ x (7) 

X 

where |0) denotes the empty state (and cr(x) means ] if +1 and { otherwise). The Hilbert 
space of the total system is the sum of the two orthogonal subspaces Ti a and TCf. Corre- 
spondingly, the Hamiltonian has the matrix structure (Vf a = XV and H a = H hop + H int ) 

H a V fa \ TT , ( A^ 



The potential term appears off-diagonal because it couples the two subspaces. The resolvent 
is also a matrix, and the resolvent equation 

1 1 1 rr 1 



E-H E-H Q E-Hq j E-H 

1 1 rr 1 



E — Hq E — Hq E — Hq 
+ 'E^Wo Vfa £^#o Vfa E^H 



reads, denoting R a = (E — H a ) 1 and Rf — (E — Hf) 1 , 

1 f R a 0\ +x f R a VR s 



E-H \ R f ) \R f VR a 
R a VR f V \ 1 



f A I a RfVRaV ) E-H' (9) 
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We want to calculate the effective propagation of the a electrons to second order in A. For 
this we need only the left upper element R e s = (^zj{) aa - It is given by 

R eS = R a + X 2 R a VR f VR cS (10) 

or in other words 

R eS ~ l = E-H a -\ 2 VR f V (11) 

Up to this point, everything was exact. Now we approximate by keeping only terms up 
to order A 2 . Moreover, since we want to see the influence of the interaction in the ground 
state only, we may calculate VRfV by its expectation with respect to /, taken in the state 
| Neel ). This gives a staggered magnetic field term 

V a+(xW 3 a(x)e iQ - x (12) 

E + e f + U V ' V ; V ' 

J x 

Finally, we want to define the effective Hamiltonian H e s by -R e ff _1 — E — H e R. To do so, we 
note that the energy E is in the a band, so \E\ < 2t « e/ + U, and thus we may replace 
E+e f +u ky 7J+U' Wftli this approximation, -R e fr _1 depends only linearly on E, and, up to 
an uninteresting shift in the energy, 

#eff = H a + V a + (x)a 3 a(x) e iQ x (13) 

€ f + U X 

is the effective Hamiltonian. Since H mt couples only the a electrons, it does not play a role in 
this argument. The coupling also produces a feedback on the / electrons, but this influences 
the a electrons only in higher order in A. This shows that a model of the type we study is 
induced, but the B field obtained by this argument is small. 



III. FORMALISM 

A. Functional integrals and mean field approximation 

To study our effective model, we use the path integral formalism [Tj|. It is equivalent 
to the Hamiltonian framework and more convenient for calculating the grand canonical 
partition function Z = Tr e - ^^ - ^ , which is thus expressed as a functional integral Z = 
over Grassmann variables ip a (x) and ip a (x) carrying a spin index a = {j, J,}. 
Configuration space is labeled by x = (t z ,x) where < r x < (3 is the usual imaginary 
time and x are vectors in the <i-dimensional cubic lattice Z d (i.e. we set the lattice constant 
equal to 1). Most of our calculations are for arbitrary d, but we are mainly interested in 
d = 2 or 3. In intermediate steps of our derivation we will implicitly assume that space 
is a finite square (cube) with a finite number L d of points, but we will eventually take the 
thermodynamic limit L — > oo. We group the spacetime argument x and the spin index a 
to a single coordinate X = (x,a), and set ip(X) = ip a (x) etc. We also use convenient short 
hand notations J2 X = Jjf dr x X! xe A and = J2 x ,a etc - 



S 



The action of the model is A = J2x $(X){d T — fi)ip(X) + H(ip, ip)dr = A2 + A4, with 
a quadratic part defining the free Green function, A2 = — ^2i XY '^{X)Gq~ 1 {X, Y)ip(Y), and 
a quartic interaction term. 

To derive mean field equations for possible superconducting states we now introduce 
auxiliary bilocal pair fields and thus decouple the fermion interaction (see e.g. 



I VA*VA exp ( 



/|A(x,r)p 



2V(X,Y) 



1 T . 



■-A(X, Y)i,{Y)i,{X) - -^(X)^(Y)A(X, Y) ) (14) 



where 



^a*^a TT dReA(X,Y)dImA(X,Y) 

VAVA= R 2*W(X,Y)\ < 15 * 

is the measure for a complex boson path integral. The product runs only over those un- 
ordered pairs {X, Y} for which V(X, Y) 7^ 0. We call the A(X, Y) Hubbard- Stratonovich 
(HS) fields. Note that the integral over A(X,Y) is convergent because (and only if) the 
pairing potential is purely attractive, i.e. V(X, Y) < whenever it is nonzero. Note also that 
due to the anticommutativity of the Grassmann variables only antisymmetric configurations, 

A(X,Y) = -A{Y,X), (16) 

contribute to the path integral. This corresponds to the Pauli principle. 
Integrating out the fermions now leads to 



VA*VAe~^ (A) (17) 



with the effective HS action 



X,Y 



where the second term come from the fermion path integral (logarithm of Pfaffian of Q 1 = 
half of logarithm of determinant of Q~ x \ and we have introduced 

g- 1 = g, 1 - s (19) 

with 

£o=(^° Jj, G (X,Y) = -G (Y,X) (20) 
the free Nambu Green function and 

S =(a* 0)' A*(X,F) = A(^X). (21) 
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In the following we write 



G F 

F* G 



(22) 



which is equal to the Nambu Green function for non-interacting fermions in an external field 
A. 

HF theory amounts to evaluating the HS path integral in Eq. (|17D using the saddle point 
method. We obtain Z ps exp[— JF(A^ F )] where A HF is the solution of the saddle point 
equation SJ r (A)/SA(X,Y) = which corresponds to the minimum of T. More explicitly 
the (complex conjugate of the) latter equation is 

A(X, Y) = V(X, Y)F(X, Y) (23) 

which together with Eq. ( [L9|) forms a self consistent system of equations and provides 
the BCS description for our model. Note that in the saddle-point approximation, 
—(^f(X)^f(Yy) = Q, so Eq. ( fL9|) can now be interpreted as the Dyson equation for our 
electron Green function Q where mean field theory gives the approximation ( |2"B| ) for the 
electron self energy S. 

The general analysis of the mean field equations is still too difficult since for general X, Y 
dependent HS field configurations A(X,Y), the evaluation of the Nambu Green function 
( p2f) is impossible in general. To proceed we therefore make the usual assumption and only 
consider HS configurations that have at least a remnant of translation invariance. Since 
for our model the free Green function contains a staggered contribution i.e. is not invariant 
under translations by one but only two sites, the simplest consistent ansatz for our HF 
equations are HS configurations A invariant by translations by two sites. 



IV. STAGGERED SUPERCONDUCTIVITY 

A. Staggered states 

We now specialize to time independent states where the superconducting order parameter 



has a uniform part and a staggered part, i.e. [20] 



A(x, y) = A^x - y) + A A (x - y)e 4Q ' y , (24) 

with Q = (jr, . . . , 7r) the antiferromagnetic vector (here and below we use spin matrix no- 
tation whenever possible). This obviously implies that Q is also a sum of a uniform and 
staggered contribution, and we can solve the Dyson equation (|TU|) by Fourier transformation. 
We use the convention 

f( x ) = e ~ ikx f( k ) > kx = u n r x - x • k 

k 

where k = (u n , k) with u n = ( 2w + 1 ^ )7r ( n £ z,) the Matsubara frequencies and k = (ki, . . . , A^) 
a momentum vector in the first Brillouin zone [— 7r, ir} d of the lattice, i.e. 
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E^E/<* 



dk = 



d d k 

(2vr) d 



Then /(Jfc) = ^e^/fa). 

With that the self-energy equation ( p3|) becomes 

A,, A (fc) = ^l/(fc-g)Vg). 



(25) 



B. T c - equations 

We first consider the region close to the critical temperature T c where the Dyson equation 
(|19|) can be linearized in the non-gauge invariant Green functions F and self energy A. 
Writing ( |i9|) as G = Go + Go * S * Go + ■ ■ ■ we obtain the linearized equation F — Go * A * Go 
i.e. 

F(x,y) = J2 G o^^i)^ l7 z 2 )Go(z 2 ,y) (26) 

21,^2 

where G (x,y) = —G (y,x) T and T means matrix transposition. With the ansatz (|24|) and 
Fourier transform this becomes (the k dependence is suppressed, F N = F N {k) etc., in the 
following) , 

Fn,a = Gon^nGon,a + Goa^%Gq AN 

+GoaA®Gon,a + Gon^aG® an (27) 

where A^(/c) = Aj^(k — Q) etc, and spin matrix multiplication is understood. Combined 
with Eq. (^5|), this gives the T c equations for our model. 

The translation-invariant and staggered parts of the free Green function, Gqn,a are 

GoN{k) = g(k)a , G oA (k) = a(k)a 3 (28a) 

where 

iuj n + fi + e(k) 



g(k) 



(iu n - E + )(iu n - E-] 



a(k) = t z— (28b) 



[we used e(k + Q) = — e(k)] and 



E±(k) = -fi± ^(k) 2 + s 2 (29) 

are the antiferromagnetic bands separated by a gap = 2s. We also set B • cr = sa 3 without 
loss of generality. 
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TABLE I. Symmetries of the order parameter. 



Symmetry 


On-site (0) 


Nearest neij 


;hbor (nn) 




Next nearest neighbor 


(nnn) 


s-wave 


% (k) = 1 


r/i(k) = cos(fci' 


+ cos(fe 2 ) 


%(k) 


= cos(fei + k%) + cos(/ci 


-Afc) 


p-wave 




%,3(k) = sin(fe! 


) ± sin(/c 2 ) 


%,7(k) 


= sin(A;i + fc 2 ) ± sin(fei 




d-wave 




774 (k) = cos(fc^ 


— cos(A? 2 ) 


%(k) 


= cos(fci + fc 2 ) — cos(fei 





To deal with the remaining dependence on the relative coordinate, x — y, we decom- 
pose the order parameter An, a in components that transform according to the irreducible 
representations of the lattice point group: 

A(x-y) = ^A, / r/ J (x-y). (30) 
J 

For a two dimensional square lattice the functions rjj we need are listed in Table 0. They 
are orthonormal, fdki]j(k)r]j/(k) = 5jj>. These functions determine the shape of the gap, 
i.e., if it is s-wave or <i-wave etc. The potential ( j3b|) has the expansion 

V(k-q) = -^^ J (k)r /j (q). (31) 
J 

where 

9i = 92 = 93 = #4 = g m 

95 = 96 = 97 = 98 = 5W (32) 

Note that the interaction V, which is a diagonal matrix in the x-basis (i.e., it acts by 
ordinary multiplication) in Eq. ( f23|) , is still diagonal in the new basis. This is because the 
coupling constants go, g nn and g nnn are the same for all on-site, for all nearest neighbor and 
for all next nearest neighbor interactions, respectively, i.e., the interaction is proportional to 
identity within each of these subspaces, and because the unitary transformation implied by 
Eq. ( p0| ) does not mix these subspaces. Then we get coupled (linear algebraic) equations for 
the 2x2 matrices (A^ t Ajj. The different channels J do not mix in the T c -equations: taking 

4 

for instance the nn interaction, we have V(k — q) = —g r/j(k)r/j(q). Inserting Eq. (j3l"]). 

j=i 

Eq. (|27|), and Eq. ( pCf ) into Eq. (pip, we get a system of equations whose coefficients contain 
integrals of the form Jc?k7/j(k)?7j/(k) ( l ) (k), where $ is a function that has the symmetries 
$(— k) = $(k) and $(/c2, k\) = $(&i,A: 2 ) because it depends on k only through e(k). The 
first symmetry implies that p-wave (J = 2,3) does not couple to s and d-wave (J = 1,4), 
and the second symmetry implies that s and d are uncoupled. 
We write 
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so that the A^j and O^j are the c-number order-parameters for the superconducting state. 
The final equations will give a critical temperature T c for each of the possible channels 
separately, and the biggest of those determines the dominating channel and the T c of the 
system. 

From now on, we suppress the indices (J, i) and write 



V = r)J, 9j = fi'Bcs 



etc.; it is clear from the above that one has to calculate the critical temperatures for all 
interesting channels and take the maximum one. This determines the dominant SC channel: 
the 'shape' r/(k), the spin structure, and the ratio of translation invariant and staggered 
parts A and of the gap. 

Every shape function 77 can be characterized by two sign factors p = +, — and £ = +, — 
defined by 

V (k-Q)=p V (k), 7/(-k)=^(k); (34) 
these determine which spin structures of the gap are compatible with the Pauli principle Eq. 



(|I~6l). Since in position space, the spatial dependence of Ajv,a(x — y) is given by the Fourier 
transform r/(x — y) = pi](x — y)e*Q'( x ~ y ) = £r/(y — x) of i](k), Eq. (|l~6|) implies the following 
conditions 



A 1A = -£A 1A , M = -& M (35a) 

and 

A 2 = -£A 3 , 2 = -£p0 3 . (35b) 

These conditions exclude many of the possible solutions of the mean field equations, e.g. all 
those where A\ 7^ with £7^—1 and 0\ 7^ with £p 7^ —1 etc. 

We now derive the T c -equations for our model. As the potential is frequency indepen- 
dent we can do the Matsubara sums analytically. After a straightforward calculation (see 
Appendix ||) we obtain 

4 M = 9bcsK p ((3)A 1a ± g BC sL p ([3) 6> M 

6>i,4 = 9bcsK p (P) 6i A ± g BCS L p (p)A 1A (36a) 

and 

^2,3 = gBcsK- p (P)A 2 ,3 =F 9bcsL-p(P)02,3 

02,3 = 9bcsK- p (P) 2 ,3 T 9bcsL- p (P)A 2!3 (36b) 

where p = ± is determined by ^(k) as above, and the functions K + (f3) etc. are given in 
Appendix ||. In fact L_ = 0, therefore half of the combinations will lead to uncoupled 
equations for the staggered and translation invariant gaps. 

In Table |I| we list the possible channels A(~x, y). We use the following notation for the 
spin matrix structure 
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TABLE II. Possible channels A(x,y). 



On-site ( p = +1 



0° s 



^(ITI)-|lt»% 
0e iQ - y (ITI)-HT))% 



Nearest neighbor ( p 



6> nn 
A™ 

{AkO) 
{AkO) 
{AkO) 



mi 

s 

mi 
P 

1111 

(/ 



0e iQ ' y {|TT) or \U)} m 

A{\m or |U)} 772,3 
0e iQ -y{|TT) or \U)} V4 

^(ITI)-HT)) + 0e lQ - y (ITI) + IIT»)^i 

A(\n)+\im + 8e iQ - y (\n)-\im)m,3 
A(\n)-\im + oe^-y(\]i) + \i]))) m 



Next nearest neighbors ( p = +1 



a nnn 

4nnn 

P 

A T n 

/Qnnn 
/Qnnn 
/Qnnn 



^(IU)-Ut))% 

^(lti)+lit»%,7 

a(\u) -\imvs 

0e iQ - y (IU)-|IT))% 
6> e iQ-y (| Tl) + UT»»76,7 

0e^(|U)-Ut))^ 
( zi + 0e iQ-y)| TT ) 7?6 

(Z\ - ©e^-y) HI) J77 



"TT> = (J °), IT1>= (° J), etc. (37) 

and ??j is short for ?7j(x — y) = J 5k e lk ^ x_y ^?7(k) i.e. the Fourier transformed functions 
defined in Table |. We use obvious symbols A° s etc. as abbreviations for the channels. 

Note that some of the gaps are a mixture of spin singlet and triplet. This comes as no 
surprise since spin rotation symmetry is broken in our model. 

We will analyze these T c -equations numerically in the next section. First, however, we 
discuss a simplified analysis which applies in the limit where the pairing across the AF gap 
can be neglected: For large s and large f3 (small temperatures), a nontrivial filling p ^ 1 is 
only possible for s ~ w ~ thus 

\K+\ » |#+| » \L+\ « M\og(Q(3) , (38) 

and are much smaller in comparison (explicit formulas for these constants are derived 

in Appendix |B]) . The constant 

M = J dk^(5(E + ) + 8(E^) V (k) 2 (39) 

has the natural interpretation as DOS N s (p) times the Fermi surface average of r/(k) 2 (shape 
of the SC gap-squared), and Q can be regarded as an effective energy cutoff of the interaction. 
Thus we get a BCS like equation 
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T c = tt e ~ 1/9BCS 



(40) 



with A = 2M, which can be trusted as long as T c is small. This already allows us to make 
predictions about the most stable channels for large s-values: the channels where A and 
are decoupled should have a negligible T c (due to the smallness of and KS) and thus 
should be irrelevant as compared to the mixed channels {AhO)''\ Moreover, for the latter 
we should always find O ~ ±zi, at least for sufficiently large s. Thus in case of dominating 
nnn attraction one should expect the SC channel (Z\&6>)p nn . In case of dominating nn 
attraction, also the shape of the gap matters, and one has to check which one leads to the 
biggest M i.e. the largest Fermi surface average of ^(k) 2 . For larger s-values, the Fermi 
surface is similar to the one of the half-filled hopping band (since the Fermi surface is given 
by E± = i.e. e(k) = ±a/p 2 — s 2 which is ~ for larger s-values at nonzero doping), and 
774 leads to the biggest M (it is easy to see why r/i cannot lead to a significant M: it is 
proportional to e and thus ~ at the Fermi surface). Thus for dominating nn attraction, 
one should expect the SC channel (ASzO)™. We confirm this result in a systematic stability 
analysis of all channels in the next section. 

The Zi&(9-structure of the dominating SC channel has a simple physical explanation: 
The density of electrons should follow the staggered magnetic field Be jC *' x : denoting as f the 
spin directions in the direction of B and J, opposite to it, the density difference (x)— P|(x) of 
spin-| and spin- J. fermions on the site x should be proportional to |B|e J< ^' x . Thus one should 
expect that stable Cooper pairs can only arise from electrons with their spins in direction of 
g e «Q x (otherwise the density of the participating electrons is very small). If both electrons 
are on the same site (on-site attraction), this is never possible and no significant SC is 
expected. In case the paired electrons are on adjacent sites (nn attraction), the preferred 
spin configuration is expected to be 1 1|) or | jf) (following the AF background), and there 
should be nearly no mixing of these two. Finally in the nnn case, the most stable Cooper 
pairs should be those with | f"f) or 1 1 j) (following the AF background). 



C. Mean field theory below T c 

The mean field theory below T c is obtained by minimizing the action, Eq. ([Tj|), for field 
configurations of the form given in Eq. (|24T). Introducing also the transformation fl3"0|) we 
can write the action as 

T = (3L d |Zlw|2 2 + j |0W ' 2 - ^Tr log (CT 1 ) . (41) 

J, a, a' 

The second term in the action can be expressed as Trlog(^ _1 ) = \(3L d log det lZ{k), 
where 1Z is the 8 x 8-matrix given by 



where R{k) = —R (—k) with 



Rn Ra \ a ( ^/ 



R -\R Q A R q n ) ' A "U^ ) (43) 
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and 

R N = [iu n - (e(k) - h)\(t q , R A = -sa 3 . (44) 

We assume that only one channel will contribute, i.e., we put all gj but one equal to zero. In 
contrast to the T c -equation, this absence of channel mixing is not enforced by a symmetry 
but is an additional assumption. Then, with Eqs. Q55D and e(k) = e(— k) = — e(k — Q), 
the determinant can be evaluated by finding the eigenvalues E c of 1Z: 

1 4 

- log det lZ(k) = - Ec(k)) (45) 

c=l 

(since every eigenvalue is 2-fold degenerate). These eigenvalues describe the electron bands. 
For the dominant channels A&0 in Table [H]> these electron bands E c are give by 

Ei t2 = ±E + E 3A = ±E_ (46a) 

with 

E± = ^e 2 + li 2 + s 2 + \A N \ 2 + |A A | 2 ± 2W 
W= v/e 2 |A A | 2 + eV + [fis + \A N \\A A \ cos(0)] 2 

e = e(k) , A N = A V {k) , A A = 6 V (k) (46b) 

where is the relative phase of A and 0, 

A = \A\e i(a+ ^, = \6\e ia (46c) 

and the bands are independent of a. [We do not write down the formulas for the other 
channels since they are irrelevant, as discussed.] 
Performing the Matsubara sum we obtain 



T = (3L d f (A, 9) (47) 



with 



/° = E lA '" l \ + q ^ ~lt f^fogil t! -*) . (48) 

a,cr' c=l 

(We note that the Matsubara sum used here was regularized: it is obtained by integrating 

\ Eu„ iun-E = l+e/^-fa ' anc ^ dropping an irrelevant infinite constant.) 
The mean field solution is now obtained by minimizing the action, 

= min/o(A 9) (49a) 

where we have to take the absolute minimum. Moreover, the chemical potential is fixed by 
filling, 
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, = -- . (49b) 

The standard mean field equations are dT jdA = dJ- ' jdO = and have, in general, several 
solutions which are not the absolute minimum and therefore have to be discarded. In the 
formulation Eq. (fl9"|), the correct solution is selected automatically. 

Clearly the minima will be invariant under global shifts of the phase a. Thus we end up 
with an simple minimization problem in the three parameters |Z\|,|(9| and (f>. We find that 
the relative phase <ft minimizes the action when W is minimized, i.e. cos0 = — sign(/x). 

With the mean field solution at hand we can now proceed to calculate physical quantities. 
The Green function is given by the matrix inverse of lZ(k). As an important example, we 
recall that the superconducting one-particle density of states is obtained as 



N(u) 



~ Im^ Jdk(G N ) aa {k,iu n -f u + . (50) 



Due to the complexity of the bands, (fR]), this expression is much more complicated than the 
standard BCS result. The differential current-voltage dependence of the tunneling current 
between a superconductor and a normal metal is related to this by 

^Uf-M^, (51 ) 



dV J v ' V duo 
where fp{E) = p $e +1 is the usual Fermi-Dirac distribution, and e is the electron charge. 



V. NUMERICAL RESULTS 

In this section we describe the results of a numerical evaluation of the mean field equations 
in the previous section. The numerical evaluation was done by performing the Brillouin zone 
integrals using a standard integration routine. Since the integrands are typically strongly 
peaked, care was taken to ensure sufficiently small discretization errors by using an adaptive 
step size method. Moreover, we solved mean field equations by minimizing the function f 
Eq. ( f48[) using standard numerical routines. 

For simplicity we neglect possible mixing of gaps assuming that the dominating gap 
transforms under an irreducible representation of the symmetry group of the lattice (note 
that below T c , this is an approximation since in general the energy dependence of the gap 
cannot be neglected as in ordinary BCS theory). We use the following parameters t = 
0.5 eV, gees = 0.05 eV, < s < 4eV motivated by high temperature superconductors 
(note that for the half filled Hubbard model, mean field theory predicts antiferromagnetic 
order with s ~ U/2). 

The results in Figs. were obtained by evaluating the T c equations (|36|) . To compare 
the stability of different gaps, we found it convenient to determine the required BCS coupling 
^bcs t° achieve a given critical temperature T c . In Fig. [j], the required gees is plotted as 
a function of the magnitude s of the AF background for all possible channels. We chose 
p = 0.9 and T c = 50 K as representative examples. One can clearly see that the mixed 
channels ASz are the only ones which are enhanced by the antiferromagnetism and always 
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dominate over the others. It is also always the same channels which dominate: (A&zO)™ 
in case of a nearest neighbor attraction, and (ASz0)p Un in the next nearest neighbor case 
[cf. Table 0|. The latter one requires a much larger coupling #bcs than the former, though. 
These results do not change qualitatively as the doping is changed over a wide parameter 
range. As s goes to zero the staggered part of the mixed channels decrease in comparison 
to the translation invariant part, and they completely decouple at s — 0. This, however, 
happens only for quite small s; for larger s they are of approximately equal magnitude. 

In Fig. the same quantity is plotted as a function of filling, p, for fixed s = 5t. One can 
clearly see the optimal doping away from half filling at which T c for fixed coupling will have 
a maximum. To illustrate the stabilizing effect of AF background, the inset in Fig. shows 
the corresponding results for s = (no AF background) for comparison. As expected, in 
this case the most stable SC occurs at half filling where the chemical potential intersects the 
van Hove singularity of the 2D DOS, but the required couplings gees to reach significant T c 
values are much larger than for nonzero s. 

We also give some typical examples of T c as a function of s, with the filling p fixed 
[Fig. 0(a)] and T c as a function of p, with s fixed [Fig. 0(b)]. We only show the dominating 
mixed channels. Obviously the antiferromagnetic background can lead to critical temper- 
atures of the same order of magnitude as experimental results for HTSC already for quite 
small coupling. Comparing our phase diagram to the experimental results for HTSC, there 
are some similarities, e.g., an optimal doping (at p « 0.9 in the nn channel for the parameters 
considered), but also some notable differences, e.g., the SC phase in Fig. 0(b) extends down 
to filling p = 1. However one should keep in mind that the doping in our model can differ 
from the exact result if for example part of the holes do not superconduct due to effects not 
included in our treatment. 

We have also tested the validity of the BCS-like formulas T c = Q exp (— l/gecsA). As 
discussed, one can expect this formula to be a good approximation for sufficiently small 
temperatures T c , but we found that it works reasonably well up to quite high T c values 
(for the parameters in Fig. 0(a) up to T c « 40 K). Q in this formula has the interpretation 
of an effective energy cut-off for the interaction, similar to the Debye frequency in the 
simple BCS model. We note that for fixed p, Q becomes independent of the shape of the 
gap and proportional to the AF band width « t 2 / s for larger s- values. For the parameters 
adequate for high temperature superconductors which we considered here, Q is huge — of 
order leV~ 10 4 K — for s = 0, but it becomes comparable to small values typical for 
standard superconductors (order 10 K) in presence of AF order. Nevertheless, without AF 
order (s = 0), one never can get a significant T c since A is too small, and the increase of T c 
with s is due to the increase of A resulting from the increase of the DOS in the AF bands: 
Since A enters T c exponentially, it overcompensates the decrease of Q with s. 

In Fig. we show the temperature dependence of the magnitude of the superconducting 
order parameters A and for parameters values realistic for HTSC [Eqs. Q49"D]. In all cases, 
\A\ is slightly bigger than \ 0\. Note that \A\ and | 0\ are nearly the same, which shows that 
the SC gap closely follows the antiferromagnetic background. This is only true for larger s- 
values. From this figure we also obtain 2|A max (T = 0)|//csT c « 4.7 (A max = max A(x, y) = 
2(|Z\| + |@|)), but we found that this parameter depends on doping and s and is not universal 
as in the BCS model. 

Another interesting quantity is the superconducting one-particle electronic density of 
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states N(u). In Fig. |5|(a) we show how N(u) changes as the temperature is lowered from 
T c , for a c? nn -wave gap [Eq. fl50l)1. Here one sees how the SC gap opens up at the Fermi 
level, which is at u = 0. The filling used in the figure is chosen to give an optimal T c for 
this value of s, and one clearly sees that the gap opens up precisely at the AF peak in N(u) 
as the temperature is lowered below T c . It is interesting that our simple model can lead 
to the rich structure in the DOS seen in the figure. In Fig. |5](b) we show the tunneling 
conductance, dl/dV, for a superconductor - normal- metal junction, which is essentially 
the same curve smeared out by positive temperature effects [Eq. (|5lQl. Here we note some 
significant features also seen in experiments pH] : (1) There is not a real gap but only a dip 
in the tunnel conductance. This is a consequence of the (i-wave shape of the gap, which 
implies that it has nodes on the Fermi surface. (2) The positions of the two peaks do not 
move much as the temperature is lowered, even though they would be expected to, because 
of the increase of the gap. It is clear, especially when compared to Fig. ^|(a), that this 
comes from thermal smearing. At low temperature, when the smearing is small, the gap is 
essentially constant as a function of temperature. Some details of the structure in the DOS 
remain visible at very low temperature. Note also that the separation of the two peaks is 
about 2|A max | (compare Figs. |]and|5|). 

The results also provide a consistency check for our assumption that the different SC 
channels do not mix below T c : for larger s values, it is always the channel A&cO™ respectively 
ASz0p nn which is dominating: for a given coupling g nn respectively g nnn , it always leads 
to a T c much bigger than the others, and also the gaps for the different channels remain 
well-separated for all temperatures. This suggests that neglecting the channel mixing was 
justified. 

VI. DISCUSSION 

In this paper we studied a simple model of superconductivity in an antiferromagnetic 
background and we analyzed in detail the two dimensional version of this model. Our 
model offers a simple account for one effect of AF correlations on superconductivity: AF 
can narrow bands and thereby boost T c . We also discussed the possible relevance of this for 
HTSC, although the connection is not firmly established. 

To summarize our results, in case of nearest neighbor (nn) attraction, we found that the 
internal structure of the gap is ci-wave, whereas for next-nearest neighbor (nnn) attraction 
it is p-wave. As soon as AF is present, the same coupling produces a much larger T c in the 
rf-wave than in the p-wave case, and on-site attraction leads only to insignificant T c values. 
The most stable SC channels correspond to paired fermions oriented parallel to the AF order 
parameter, and in case the paired fermions are on adjacent sites (nn) they have opposite 
spins and correspond to a mixture of spin singlet and triplet. The latter is possible because 
spin rotation symmetry is broken. In the <i-wave channel, we found that for reasonably 
weak attraction, switching on the AF order increased T c to the right order of magnitude for 
HTSC. 

We stress that all the results presented in this paper are within a mean field approxi- 
mation. We expect fluctuation corrections to the curves in all the figures, and the results 
should thus be regarded only as order of magnitude estimates demonstrating the qualitative 
features of the model. 



19 



A comment about the apparent long range order of our mean field solution is in order. 
The mean field analysis done here allows to determine only the magnitude of the SC order 
parameter but not its phase, and our calculated T c rather corresponds to the mean-field 
transition temperature T c0 . Thus nontrivial solutions of our mean field equations do not 
contradict the nonexistence of continuous symmetry breaking in two dimensions. A complete 
analysis of SC in two dimensions also has to take into account phase fluctuations |22[ and 
is beyond the scope of this paper. 

To conclude, we presented a simple model for electrons in an AF background, and studied 
superconductivity in this model. We discussed how this model may describe certain aspects 
of the physics of antiferromagnetically correlated fermions, as present e.g. in HTSC. Some 
of our results are qualitatively similar to certain basic properties of HTSC. The methods 
presented here should be useful also for other similar and more complicated models, in 
particular for cases where several effective bands contribute to superconductivity. 
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APPENDIX A: ON STAGGERED STATES 

In this Appendix we show how to evaluate the inverse for Green functions etc. describing 
staggered states. Note that the quantities A(x, y) etc. below are matrices and the matrix 
product is understood. 

To get the inverse of A(x, y) = A N (x — y) + Aa(x — y)e tQy , we solve the equation 



^2(A N (x -z) + A A (x - z)e i 



x(B N (z -y) + B A {z - y)e iQy ) = S x>y , (Al) 
which, after Fourier transformation, can be written in a matrix form 

A N A A \ (B N B A \ (1 0\ 

Ai A%){Bj B%) = {» l) < A2 > 
where A^ (k) = A N (k — Q), etc. The solutions are (suppressing the argument k) 

B N = (A N - ^(Ag)- 1 ^)- 1 

B A = {A Q A -A Q N A' A 1 A N y 1 (A3) 
If all of An, A a , Bn, Ba commute, 

A Q 

B = N 
N A N A Q N - A A A Q A 

B A = ^ q (A4) 

AnA n — AaA a 

In particular, this gives the free Green function ( p8| ) as the inverse of R = Gq 1 has the 
translation invariant and staggered parts Rn,a given in Eq. 
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APPENDIX B: T c EQUATIONS 



Here we give some details of our derivation of the T c equation in Section [IV B . 
We use the following notation for functions of k: 

fQ(k) = f(k-Q), /(*) = -/(-*) 

We also write Gjv(^) = g(k)&Q, Ga(A;) = a(k)cr s and use that rj Q {k) = rj(k — Q) = pr)(k), 
where p = ±1 depending on the parity of A(k). 
With that we obtain the equations 

^1,4 = -£bcs ^2({ggV 2 + ad Q rjr] Q )A 1A 

k 

±{ga Q if + agr]r] Q )6 1A 
0i,4 = -£bcs ((^V + adr]r] Q ) O x 

k 

±(gdr]r] + ag® r\rp") A\,\ 

and 

^2,3 = -fi-Bcs ^((ggv 2 - aa Q r]7] Q )A 2) z 

k 

T(ga Q r] 2 - agr]7] Q ) 6> 2 , 3 ) 
02,3 = -^bcs ^2 ((99 Qf l 2 ~ adr]i] Q ) 6> 2 , 3 

k 

T(gdr] 2 - ag Q r]ri Q )A 2t3 

Doing the Matsubara sums and defining 

v(n Jy 1 tanh(/?E/2) 

W E ) = pl^^&- 2E (B1) 

results in Eq. fl56|) with the constants 



s 



2 



-( X p{E + ) - Xp{E-)Wj 



2fiw 

K + (f3) = jdk (±(xfi(E + ) + Xp{E-M 



e 



2 



( X p(E + ) - Xfs(E-))r]j 



2fiw 

21 



111 



L_(/3) = (B2) 



where w = \J e 2 + s 2 and E± = —fx ± w. Note that all these expressions have well-defined 
limits fi —>■ 0. 

It is shown below (and can be checked numerically) that the f3 dependence of these 
integrals for large f3 is in a very good approximation given by (these formulas becomes exact 
for P — ► oo ) 



K + (P) = M\og(Q^P), K_{(3) 



a 



s 2 



K + (p) = — Mlog(fi (fe) /3) + a, K„{(3) = a + a (B3) 
L + (/3) = -Mlog(^ c >/3) 

with M Eq. ( |55| ) and a = 073 and 5 = are, in a very good approximation, independent of 
(3 for large (5 [the latter can be checked numerically; the energy scales il^ a ' b ' c ~> can be obtained 
by calculating the integrals K + (f3) etc. for one (large) /3-value numerically and comparing 
with (|B3|).l This yields the T c -equations of the form Eq. (f40|). 

The step going from (|B2|) to (|B3l ) uses that functions F(E) reasonably smooth close to 
E — 0, the /3-dependence of 

HF.^JiEF^^^ (B4, 
for large /? is approximated well by 

J(F,/3) = F(0)bg(fi F /3) (B5) 

To see this, define 

then J(F, /?) = F(0) log(fi(F, /?)/?). Now, with x = /3E/2, 

dn n / r dx F(2 X /f3) 1 



d/3 p\J F(0) 2cosh 2 (a;) 

Taylor expanding F, we see that dQ/dp = 0(1/ P 3 ) so we can replace fl(F, P) by Qp 
lim Q(F, P) up to terms of order 1//3 2 . 
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FIG. 1. Comparison of stability of different SC channels in 2D model: required BCS coupling 
ffBCS to achieve a given superconducting critical temperature T c = 100K(t/leV), as a function of 
the magnitude of the antiferromagnetic background s at fixed filling p for different shapes of the 
gap, which are listed in Table |. (For a hopping parameter t = 0.5 eV adequate for HTSC this 
corresponds to T c = 50 K). The filling p = 0.9 considered here is representative, (a): A&i 0-channel 
i.e. mixture of staggered and translation invariant gaps which we find to dominate throughout, (b): 
Z\-channel i.e. pure translation invariant gap. (c): (9-channel, i.e. pure staggered gap. Note the 
different scales of the gees-axes in the different plots. There is also an s nn -channel (not shown in 
the figure) in both the mixed and purely staggered channel which requires an even higher coupling. 
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FIG. 2. Required coupling <?bcs to reach a critical temperature T c = 100K(t/leV) as a function 
of the filling p for fixed magnitude of the AF background, s = 5i. The different curves correspond 
to different shapes of the gap, for the mixed channels, Z\&(9. The purely translation invariant and 
staggered solutions require a much stronger coupling. The inset shows the corresponding result for 
s = (no AF). 
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FIG. 3. Critical temperature T c as a function of magnitude of antiferromagnetic background 
s and filling p, for dominating SC channel A&O and different shapes of the gap. Used parameters 
here are: t = 0.5 eV and gscs = 0.05 eV. (a) p = 0.9 fixed (b) s = leV fixed. Labeling of results 
for different shapes of the gap as in Fig. |]. 
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FIG. 4. Magnitude of the superconducting gaps, \A\ and |(9|, as a function of temperature 
T for dominating channels A&0. Two neighboring curves show the translation invariant (A) and 
staggered (0) parts of the gap. In all cases, \A\ is bigger than | 0\. Note that \A\ and | 0\ are nearly 
the same here, which shows that the SC gap closely follows the antiferromagnetic background, as 
discussed in the text (they become identical for s — > oo). The parameters are p = 0.9, s = leV, 
t = 0.5 eV and g = 0.05 eV. 
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FIG. 5. (a) Superconducting density of states for three different temperatures T = 75 K, 65 K 
and IK. The parameters are the same as in Fig. ||, with a T c ~ 75 K. The chemical potential is 
at <jj = 0. (b) Tunneling conductance dl/dV for a superconductor - normal-metal junction as a 
function of the applied voltage for a sequence of temperatures T = 75, 65, 55 . . . , 5 K and 1 K. Note 
that most of the fine structures in the density of states is washed out by temperature effects. 
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